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ABSTRACT 



Context. The nature of the several microlensing events observed by the MACHO team towards the LMC still remains controversial. 
Low-mass substellar objects and stars with masses larger than ~ 1 Mq have been ruled out as major components of a MACHO 
galactic halo, while stars of half solar masses are the most probable candidates. 

Aims. In this paper we assess jointly the relative contributions of both red dwarfs and white dwarfs to the mass budget of the galactic 
halo. 

Methods. In doing so we use a Monte Carlo simulator which incorporates up-to-date evolutionary sequences of both red dwarfs 
and white dwarfs as well as detailed descriptions of our Galaxy and of the LMC. We explore the complete mass range between 0.08 
and 1 Mq as possible microlensing candidates and we compare the synthetic populations obtained with our simulator with the results 
obtained by the MACHO and EROS experiments. 

Results. Our results indicate that, despite that the contribution of the red dwarf population increases by a factor of 2 the value of 
the optical depth obtained when taking into account the white dwarf population alone, it is still insufficient to explain the number of 
events claimed by the MACHO team. 

Conclusions. Finally, we find that the contribution to the halo dark matter of the whole population under study is smaller than 10% at 
the 95% conficence level. 

Key words, stars: red dwarfs — stars: white dwarfs — stars: luminosity function, mass function — Galaxy: stellar content — Galaxy: 
stellar content — Galaxy: structure — Galaxy: halo 



1. Introduction 

Several cosmological observations show compelling evidence 
that baryons represent only a small fraction of the total matter 
in our Universe and that non-baryonic dark matter dominates 
over baryons. Within the standard cosmological model this nat- 
urally leads to the existence of a sort of unknown energy, the 
dark energy, which is the predominant component, Qa ^ 0.72, 
along with the dark matter component, Qu - 0.27, whereas the 
baryonic component is just Ob - 0.044. Moreover, most of the 
baryons are non-luminous, since Q,^ ^ 0.005. For the case of 
our own Galaxy it has been found that the virial mass out to 100 
kpc is M « 10'^ Mq while the baryonic mass in the form of stars 
is w 7 X 10'" Mq, which means that for the Milky Way, the 
baryon fraction is at most 8% (Klypin et al. 2007). This prob- 
lem is known as the missing baryon problem — see the review 
of Silk (2007) for a complete, interesting and recent discussion 
of this issue — and it is critical in our understanding of how the 
Galaxy (an by extension other galaxies) were formed and will 
ultimately evolve. To solve this problem, three alternatives can 
be envisaged: either these baryons are in the outer regions of our 
Galaxy, or, perhaps, they never were present in the protogalaxy 
or, finally, they may have been ejected from the Milky Way. The 
most promising explanation and the currently favored one is the 
first of these options. 
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Since the pioneering proposal of Paczyriski (1986) that grav- 
itational microlensing could be used to clarify the nature of 
galactic dark matter, considerable observational and theoreti- 
cal efforts have invested in this issue. Among the most likely 
candidates for building up the baryonic dark matter density 
are massive baryonic halo objects, or MACHOs. From the the- 
oretical point of view it has been suggested that MACHOs 
could be planets (M ~ 10"^ Mq), brown and red dwarfs (with 
masses ranging from ~ 0.01 to ~ 1 Mq), primordial black holes 
(M ^ 10"'^ Mq), molecular clumps (M ~ 1 Mq) and old white 
dwarfs (M ~ 0.6 Mq). From the observational perspective the 
MACHO (Alcocket al. 1997, 2000), EROS (Lasserre et al. 2001; 
Goldman et al. 2002; Tisserand et al. 2007), OGLE (Udalski et 
al. 1994), MOA (Muraki et al. 1999) and SuperMACHO (Becker 
et al. 2005) teams have monitored millions of stars during several 
years in both the Large Magellanic Cloud (LMC) and the Small 
Magellanic Cloud (SMC) to search for microlensing events. 
However, only the MACHO team and the EROS group have 
reported their results. Despite that initially some divergences 
appeared between their results, agreement has been reached in 
some important points. For instance, no microlensing candidates 
have been found by the MACHO team or the EROS group with 
event durations between a few hours and 20 days. This means 
that the Galactic halo could contain no more than a 10% of dark 
objects in the mass range 10"^ < M/Mq < 10^^. This rules 
out planets and brown dwarfs as the major contributors to the 
mass budget of the dark halo. Moreover the MACHO collabora- 
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tion has succeeded in revealing ~ 15 microlensing events during 
their 5.7 yr analysis of 11.9 million stars in the LMC (Alcock 
et al. 2000). In their analysis they derived an optical depth to- 
wards the LMC of T = 1-2^Q 3 x 10"^ or, equivalently, an halo 
fraction 0.08 < / < 0.50 at the 95% confidence level with a 
MACHO mass in the range 0.15 < M/M© < 0.50, depending 
on the halo model. On the other hand, the negative results found 
by the EROS collaboration provide an upper limit. Specifically, 
the EROS team has presented an analysis of a subsample of 
bright stars belonging to the LMC (Tisserand et al. 2007), in 
order to minimize the source confunsion and blending effects. 
Their results imply that the optical depth towards the LMC is 
T < 0.36 X 10"' at the 95% confidence level, corresponding to 
a fraction of halo mass of less than 7%. This result is 4 times 
smaller than that obtained by the MACHO team. Consequently, 
further discussion has been opened concerning the location and 
nature of the lenses. In particular, recent LMC models have been 
used in order to ascertain if possible asymmetries in the space 
distribution of the microlensed stars are assimilable to Galactic 
halo objects or LMC ones. In fact, different studies show that 
a sizeable fraction of the microlensed events are due to LMC 
self-lensing (Sahu 1994; Gyuk, Dalai & Griest 2000; Calchi 
Novati et al. 2006). Moreover, a full variety of possible explana- 
tions have been proposed to reproduce the microlensing events: 
tidal debris or a dwarf galaxy toward the LMC (Zhao 1998), a 
Galactic extended shroud population of white dwarfs (Gates & 
Gyuk 2001), blending effects (Belokurov, Evans & Le Du 2003, 
2004), spatially varying mass functions (Kerins & Evans 1998, 
Rahvar 2005), and other explanations (Holopainen et al. 2006). 
However, all of these proposals have been received with some 
criticism because none of them fully explains the observed mi- 
crolensing results. 

Other observational pieces of evidence are added to conform 
the present puzzle, such as the search for very faint objects in 
the Hubble Deep Field or the search for the microlensing events 
towards the Galactic bulge or towards very crowded fields like 
M3 1 . Particularly, the Hubble Deep Field has provided us with 
an opportunity to test the contribution of white dwarfs to the 
Galactic dark matter. Ibata et al. (1999) and, most recently, Kilic 
et al. (2005) have claimed the detection of some white dwarf 
candidates among several faint blue objects. These objects ex- 
hibit significant proper motion and, thus, are assumed to belong 
to the thick disk or the halo populations. Despite the increasing 
number of surveys searching for white dwarfs — like the Sloan 
Digital Sky Survey (Eisenstein et al. 2006), the 2 Micron All 
Sky Survey (Cutri et al. 2003), the SuperCosmos Sky Survey 
(Hambly, Irwin & MacGillivray 2001), the 2dF QSO Redshift 
Survey (Vermes et al. 2002), and other observational searches 
(Knox et al. 1999; Ibata et al. 1999; Oppenheimer et al. 2001; 
Majewski & Siegel 2002; Nelson et al. 2002) — and the numer- 
ous theoretical studies (Isern et al. 1998; Reyle et al. 200 1 ; Flynn 
et al. 2003; Hansen & Liebert 2003; Garci'a-Berro et al. 2004), 
the problem still remains open. In addition, despite the fact that, 
for obvious reasons, no information about the halo could be de- 
rived from the microlensing events towards the Galactic bulge, 
the results obtained so far indicate the primordial role played by 
low-mass stars rather than other objects. 

In two previous papers we have extensively analyzed the 
role played by the carbon-oxygen (CO) white dwarf population 
(Garcfa-Berro et al. 2004) as well as the contribution of oxygen- 
neon (ONe) white dwarfs (Camacho et al. 2007). In these pa- 
pers we have performed a thorough study of a wide range of 
Galactic inputs, including different initial mass functions and 
halo ages, and several density profiles corresponding to differ- 



ent halo models. The calculations reported in these two papers 
have shown that a sizeable fraction of the halo dark matter can- 
not be locked in the form of white dwarfs. Specifically, we have 
found that this contribution is of the order of a modest 5% in 
the most optimistic case, that it is mainly due to old CO white 
dwarfs with hydrogen-rich atmospheres and that the contribu- 
tion of ONe white dwarfs is minor, because although ONe white 
dwarfs can reach faint magnitudes faster than CO white dwarfs, 
their contribution is heavily suppressed by the initial mass func- 
tion. 

In this paper we analyze in a comprehensive way a wide 
range of masses 0.08 < M/Mq < 10 susceptible to produce 
microlensing events towards the LMC and, thus, to contribute to 
the halo dark matter. The full range of masses studied here rep- 
resents nearly the 90% of the stellar content, including the red 
dwarf regime (M > 0.075 Mq) the CO white dwarf population 
and the population of massive ONe white dwarfs. The paper is 
organized as follows. In Sect. 2 we summarize the main ingredi- 
ents of our Monte Carlo code and other basic assumptions and 
procedures necessary to evaluate the microlensing optical depth 
towards the LMC. Section 3 is devoted to describe our main re- 
sults, including the contribution of red dwarfs and white dwarfs 
to the microlensing optical depth towards the LMC, and we com- 
pare our results to those of the MACHO and EROS teams. In this 
section we also estimate the probability that a microlensing event 
could be assigned to one or another of the populations under 
study and we discuss the contribution of red and white dwarfs to 
the baryonic content of the Galaxy. Finally, in Sect. 4 our major 
findings are summarized and we draw our conclusions. 

2. The model 

2.1. Building the sample 

A detailed description of our Monte Carlo simulator has been 
already presented in Torres et al. (2002), Garcfa-Berro et al. 
(2004) and Camacho et al. (2007). Therefore, we will only sum- 
marize here the most important inputs. The basic ingredient of 
any Monte Carlo simulator is a random number generator algo- 
rithm which must ensure a non-correlated sequence and a good 
set of statistical properties. We have used a random number gen- 
erator algorithm (James 1990) which provides a uniform prob- 
ability density within the interval (0, 1) and ensures a repetion 
period of ^ 10^^, which is enough for our purposes. Each one 
of the Monte Carlo simulations presented here consists of an en- 
semble of ~ 5 X 10'* independent realizations of the synthetic star 
population, for which the average of any observational quantity 
along with its corresponding standard deviation were computed. 
Here the standard deviation means the ensemble mean of the 
sample dispersions for a typical sample. 

The galactic halo has been modelled assuming a spherically 
symmetric halo. In particular the model used here is the typi- 
cal isothermal sphere of radius 5 kpc also called the "S-model", 
which has been extensively used by the MACHO collaboration 
(Alcock et al. 2000; Griest 1991). The position of each syn- 
thetic star is randomly chosen according to this density profile. 
We have not used other profiles, such as the exponential power- 
law model, the Navarro, French & White (1997) density pro- 
file and others because their inclusion in the analysis presented 
here does not provide significant variations in the final results 
(Garcfa-Berro et al. 2004). Likewise, this reasoning is extensi- 
ble to the non-standard galactic halo models, such as flattened 
profiles, oblate halo models, and others, which are beyond the 
scope of the present study. 
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The mass distribution of synthetic stars has been computed 
according to two different initial mass functions, the standard 
initial mass function of Scalo (1998) and the biased log-normal 
initial mass function proposed by Adams & Laughlin (1996), 
the later being representative of other non-conventional initial 
mass functions such as the one of Chabrier et al. (1996). Note, 
however, that these biased initial mass functions seem to be in- 
compatible with the observed properties of the halo white dwarf 
population (Isern et al. 1998; Garcfa-Berro et al. 2004), with the 
contribution of thermonuclear supernovae to the metallicity of 
the Galactic halo (Canal et al. 1997), and with the observations 
of galactic halos in deep galaxy surveys (Chariot & Silk 1995). 
Nevertheless, for the sake of completeness, we prefer to include 
in the present analysis a representative example of these biased 
mass function in order to illustrate the role played by the red 
dwarf population when a biased initial mass function is adopted. 

The main sequence mass is obtained by drawing a pseudo- 
random number according to the adopted initial mass function. 
The time at which a star was born is randomly choosen, pre- 
viously assuming that the halo was formed 14 Gyr ago in an 
intense burst of star formation of duration ~ 1 Gyr. The main- 
sequence lifetime as a function of the mass in the main squence 
is that of Iben & LaughUn (1989). Once the mass of a syn- 
thetic star is chosen, the main-sequence lifetime is obtained 
and from it we know which stars are able to evolve to white 
dwarfs or which ones remain in the main sequence as red dwarfs. 
We have considered red dwarfs to have masses in the range 
0.08 < M/Mq < 1. For these stars we have adopted the evo- 
lutionary models of Baraffe et al. (1998). Stars with such small 
masses have very large main-sequence lifetimes. Consequently, 
no post-main-sequence evolutionary tracks of these stars were 
needed. For those stars which have had time enough to enter 
in the white dwarf cooling track, and given a set of theoretical 
cooling sequences and the initial to final mass relationship (Iben 
& Laughlin 1989) their luminosities, effective temperatures and 
colors were obtained. The cooling sequences adopted here de- 
pend on the mass of the white dwarf. White dwarfs with masses 
smaller than Mwd = 1.1 Mq are expected to have a CO core and, 
consequently, we adopt for them the cooUng tracks of Salaris et 
al. (2000). White dwarfs with masses larger than Mwd = 1-1 Mq 
most probably have ONe cores and for these white dwarfs we 
adopt the most recent cooling sequences of Althaus et al. (2007). 
Both sets of coohng sequences incorporate the most accurate 
physical inputs for the stellar interior (including neutrinos, crys- 
tallization, phase separation and Debye cooling) and reproduce 
the blue turn at low luminosities (Hansen 1998). 

The kinematical properties of the halo population have been 
modeled according to gaussian laws (Binney & Tremaine 1987) 
with radial and tangential velocity dispersions accordingly re- 
lated by the Jean's equation and fulfilling the flat rotation curve 
of our Galaxy. We have adopted standard values for the circu- 
lar velocity Vc = 220 km/s as well as for the peculiar velocity 
of the Sun (Uq, Vq, Wq) = (10.0, 15.0, 8.0) km/s (Dehnen & 
Binney 1998). Additionally, we have discarded those stars with 
velocities smaller than 250 km/s given that they are not con- 
sidered as halo members. Besides we have rejected stars with 
velocities larger than 750 km/s, because they would have veloci- 
ties exceeding 1 .5 times the escape velocity. Finally, since white 
dwarfs usually do not have determinations of the radial com- 
ponent of the velocity, the radial velocity is eliminated when a 
comparison with the observational data is needed. 

Finally, and in order to compare the simulated results with 
the observational ones, a normalization criterion should be used. 
We have proceeded as in our previous papers (Camacho et al. 



2007; Garcfa-Berro et al. 2004). That is, we have normalized our 
simulations to the local density of halo white dwarfs obtained 
from the halo white dwarf luminosity function of Torres et al. 
(1998) taking into account the new halo white dwarf candidates 
in the SDSS Stripe 82 (Vidrih et al. 2007). 

2.2. The LMC microlensing model 

In order to mimic the microlensing experiments towards the 
LMC we have simulated it following closely the detailed LMC 
descriptions of Gyuk et al. (2000) and Kallivayalil et al. (2006). 
Our model takes into account, among other parameters, the scale 
length and scale height of the LMC, its inclination and its kine- 
matical properties. This model provides us with a synthetic pop- 
ulation of stars representative of the monitored point sources. 
Afterwards we have evaluated which star of the galactic halo 
could be responsible of a microlensing event. We have only con- 
sidered stars fulfilling a series of conditions. First of all the lens- 
ing star should be fainter than a certain magnitude limit. In a 
second step we have checked if the lens is inside the Einstein 
tube of the monitored star. That is, if the angular distance be- 
tween the lens and monitored star is smaller than the Einstein 
radius. We recall here that the Einstein radius is given by 
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where Dqs is the observer-source distance and x = Dql/Dos- 
Finally, we filter those stars which are candidates to produce a 
microlensing event with the detection efficiency function, e(f,), 
where f, is the Einstein ring diameter crossing time. The detec- 
tion efficiency depends on the particular characteristics of the 
experiment. In our case we have reproduced the MACHO and 
EROS experiments. Specifically, for the MACHO collaboration 
we have taken 1.1 x 10^ stars during 5.7 yr and over 13.4 deg^, 
whereas the detection efficiency has been modelled according to: 



sit) 



0.43e-W'/r..))"Vo.87^ f >r„ 

43e-|ln{f/7-„)|"Vll.l6 f<T 



(2) 



where = 250 days. This expression provides a good fit to 
the results of Alcock et al. (2000). For the EROS experiment we 
have used 0.7 x 10^ stars over a wider field of 84 deg^ and over 
a period of 6.7 yr. Regarding the detection efficiency we have 
adopted a numerical fit to the results of Tisserand et al. (2007). 

For all the simulations presented here we extract the param- 
eters relevant to characterize the microlensing experiments. A 
complete description of the various parameters which have im- 
portance in discussing gravitational microlensing can be found 
in MoUerach & Roulet (2002) and Schneider et al. (2004). 
Among these parameters perhaps the most important one for our 
purposes is the optical depth, r, which measures the probability 
that at a given time a star is magnified by a lens by more than a 
factor of 1.34. From an observational point of view an estimate 
of this parameter can be obtained using the expression (Alcock 
et al. 2000): 
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where E is the total exposure in star-years. The optical depth is 
independent of the lens motion and mass distribution. However, 
since the experiments measure the number of events and their 
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durations, additional information can be obtained using the mi- 
crolensing rate, F, and its distribution as function of the event 
durations. This parameter represents nothing else but the flux of 
lenses inside the microlensing tube. Finally, an estimate of the 
expected number of events can be done using the expression 
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3. Results 

3.1. Optical depth towards the LMC 

As previously mentioned, the optical depth provides us with the 
most immediate and simplest information about the microlens- 
ing experiments. Even so, the optical depth plays a key point in 
our analysis since it provides us with rich information about the 
Galactic halo and the role of the dark matter. We will compare 
here the results obtained using our Monte Carlo simulator with 
those obtained by the MACHO collaboration. In Fig. 1 we show 
the contribution to the optical depth of the diff"erent simulated 
populations for the two initial mass functions under study as a 
function of the adopted magnitude cut. Our simulations have 
been normalized to the value derived by Alcock et al. (2000), 
To - 1.2 X 10"^. The white dwarf populations are represented 
by solid and open squares for the CO and ONe white dwarfs, re- 
spectively. Open triangles show the contribution of red dwarfs, 
while the contribution of the three populations is represented by 
open circles. The first remarkable result is that for the standard 
initial mass function (top panel of Fig. 1) the combined contri- 
bution of red dwarfs and white dwarfs is at most one third of 
the observed optical depth and this is obtained when a totally 
unrealistic magnitude cut is adopted, my ~ 15™^^. Note as well 
that there is a clearly decreasing trend as the adopted magnitude 
cut increases. This can be easily understood, since fewer objects 
contribute to the optical depth as the magnitude cut increases. 
Furthermore, the slope is different for the three types of objects 
considered here. For instance, the contribution of red dwarfs de- 
creases faster than the contribution of CO white dwarfs as the 
magnitude cut increases which, in turn, decreases faster than the 
contribution of ONe white dwarfs. This reflects the fact that, in 
general, red dwarfs are bighter than regular CO white dwarfs. 
Finally, ONe white dwarfs cool very fast (Althaus et al. 2007) 
and thus, for realistic halo ages, are faint objects. Consequently, 
their contribution remains almost constant. It is interesting to 
note as well that the optical depth obtained for the whole popu- 
lation almost doubles that obtained when only the white dwarf 
population is taken into account. However, in agreement with 
other studies, the value obtained here still remains far below the 
observed one and it is thus clear that the Galactic halo population 
is not enough to explain the results of the MACHO team. Even 
the alternative initial mass functions, of which that of Adams & 
Laughlin (1996) is a representative example, predict that (t/tq) 
is, in the best of the cases, 0.3 and this happens for unrealistic 
magnitude cuts — see the bottom panel of Fig. 1 . The results ob- 
tained using the non-standard initial mass function deserve two 
additional comments. In particular, it is interesting to note that 
the largest contribution in this case is that of the CO white dwarf 
population and, secondly, that the contribution of red dwarfs is 
totally negligible for magnitude cuts larger than my ~ 23™"^, 
which is a reasonable value for current surveys. Both facts are 
natural since this initial mass function has been tailored to pro- 
duce a large population of 0.5 Mq white dwarfs. In summary, 
the fraction of the microlensing optical depth which a stellar halo 
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Fig. 1. Microlensing optical depth towards the LMC as a func- 
tion of the limiting magnitude. Solid and open squares repre- 
sent the CO and ONe white dwarf populations, respectively. Red 
dwarfs are represented using open triangles, while the whole 
population is shown using open circles. 



can account for in the range of masses under study here is at most 
one third of that observed by the MACHO team, independently 
of the adopted initial mass function. 

In Table 1 we have summarized the average values for the 
different parameters of the whole population for the two ini- 
tial mass functions under study and for four magnitude cuts. 
Speciflcally, in Table 1 we present the number of microlens- 
ing events, the average mass of the microlenses, their average 
proper motion, distance and tangential velocity, the correspond- 
ing Einstein crossing times and, finally, the contribution to the 
microlensing optical depth. It is clear that some of these pa- 
rameters are to some extent dependent of the magnitude cut. 
For instance, the average distance of the sample increases as the 
magnitude cut is larger. This is a natural consequence of select- 
ing more distant objects which, in turn, impUes longer Einstein 
crossing times. This behavior is independent of the assumed ini- 
tial mass function. However, there are characteristics which do 
not change as the magnitude cut increases. For instance, this is 
the case of the expected number of events or the average mass 
of the sample. In the case of a standard initial mass function 
no more than one microlensing event should be expected at the 
Icr confidence level, while for a log-normal mass function up 
to 5 events might be expected. In any case the expected number 
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Standard AL 



Magnitude 


17.5 


22.5 


27.5 


32.5 


17.5 


22.5 


27.5 


32.5 


<A^WD> 


± 1 


0+1 


± 1 


± 1 


3 ± 2 


2 ± 1 


1 + 1 


0+1 


(m) (M/Mq) 


0.421 


0.411 


0.427 


0.443 


0.638 


0.636 


0.640 


0.684 


(M) ("yr-') 


0.020 


0.017 


0.008 


0.004 


0.036 


0.025 


0.011 


0.003 


(d) (kpc) 


2.48 


3.79 


6.62 


13.08 


1.39 


2.15 


5.13 


19.6 


<V,an>(kms-') 


240 


247 


262 


241 


240 


252 


263 


261 


<fE> (d) 


41.2 


49.3 


63.3 


82.8 


34.7 


46.6 


76.4 


126.8 


<T/To> 


0.283 


0.214 


0.139 


0.055 


0.302 


0.204 


0.140 


0.129 



Table 1. Summary of the results obtained for the whole simulated population of microlenses towards the LMC for an age of the 
halo of 14 Gyr, different model initial mass functions, and several magnitude cuts. 



of microlensing events is far from the 17 events claimed by the 
MACHO experiment. 

As already mentioned, the average mass of the microlenses 
depends on the initial mass function. To investigate this further 
in Fig. 2 the contribution to the optical depth as a function of the 
mass of the lens object for both initial mass functions is shown. 
The results obtained with our Monte Carlo simulator clearly 
show that for a standard initial mass function — top panel of 
Fig. 2 — there are two peaks centered at masses ~ 0.3 Mq and 
~ 0.6 Mq, respectively. These masses correspond to the average 
masses of red dwarfs and CO white dwarfs, respectively. Also 
noticeable is that the contribution of ONe white dwarfs can only 
be appreciated as an extended tail in the case of a standard initial 
mass function. This is in full agreement with recent studies about 
the distribution of masses of the white dwarf population (Finley 
et al. 1997; Liebert et al. 2005), which show the existence of a 
narrow sharp peak near 0.6 Mq, with a tail extending towards 
larger masses, with several white dwarfs with spectroscopically 
determined masses within the interval comprised between 1.0 
and 1.2 Mq. The situation is different for the non-standard ini- 
tial mass function, which is shown in the bottom panel of Fig. 
2. The log-normal initial mass function considered here cannot 
produce red dwarfs with masses below ~ 0.45 Mq and thus the 
peak at 0.3 Mq previously found is absent in this case. 

A more detailed analysis of the role played by the red dwarf 
population can be done. To this end in Table 2 we have summa- 
rized the average values for the different parameters of the pop- 
ulation of red dwarfs for both initial mass functions. Similar sets 
of data for the white dwarf population can be found in our pre- 
vious studies (Camacho et al. 2007, Garcia-Berro et al. 2004). 
As can be seen from Table 2, the red dwarf population roughly 
contributes a 10% to the observed MACHO optical depth for 
a standard initial mass function. It is also important to discuss 
the other parameters shown in Table 2. For instance, the average 
mass clearly decreases when the magnitude cut increases, which 
is the opposite to what occurs for the white dwarf population. 
This can be easily explained by taking into account that the more 
massive the red dwarf is the brighter is and, hence, we should ex- 
pect less massive objects for larger magnitude cuts. This result 
is reinforced by the fact that the average distance increases for 
increasing magnitude cuts. Moreover, since the average tangen- 
tial velocity remains constant, the combined effect of an average 
decreasing mass and an average increasing distance is that the 
Einstein crossing time remains practically constant. The charac- 
teristics of the red dwarf population are dramatically different 
when the initial mass function of Adams & Laughlin (1996) is 
used. In this case the production of low-mass red dwarfs is heav- 
ily suppressed. Accordingly, in our simulations we have not pro- 
duced red dwarfs with masses smaller than ~ 0.45 Mq. Thus, 
since in average the masses are larger we also find brighter stars. 
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Fig. 2. Contribution to the optical depth as a function of the lens 
mass. 



Hence, we expect no contribution at all for magnitude cuts above 
my - 22.5™^^, while for brighter magnitude cuts the average 
mass expected is ~ 0.7 Mq, which is even larger than the ex- 
pected value for CO white dwarfs. 

As already seen, for the standard initial mass function, a 
double-peaked profile is found, but the peak amplitude deserves 
a detailed analysis. The ratio of the contribution to the optical 
depth of a typical red dwarf with respect to the contribution of a 
typical CO white dwarf is 



TCO 



fRD s(ico) MrdD^^ s(ico) 



tco e(?RD) ~ \ McoD^^ eiijiD) 



(5) 
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Table 2. Summary of the results obtained for the population of red dwarf microlenses towards the LMC for an age of the halo of 
14 Gyr, different model initial mass functions, and several magnitude cuts. 
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Fig. 3. Fraction of microlenses with respect to the whole popu- 
lation, as a function of the magnitude cut. 



This value depends on the adopted initial mass function and on 
the magnitude cut used. For the standard initial mass function 
and a magnitude cut of 22.5"^^ the average values of the mass 
and distance of a red dwarf are, respectively, ~ 0.228 Mq and 
4.88 kpc, while for a typical CO white dwarf the average mass is 
~ 0.568 Mq and the average distance is 3. 14 kpc. With these val- 
ues the ratio of the optical depths turns out to be trd/t"co ~ 0.9. 
Thus, although the sample space of red dwarfs and white dwarfs 
is different the contribution to the optical depth per object re- 
sults to be the same. For instance, the red dwarf population has 
a smaller average mass than white dwarfs but, conversely, the 
average distance is larger. Therefore, the total contribution of 



these two populations results to be merely an account of the 
number of microlenses. In order to clarify this issue we have 
evaluated the fraction of microlenses due to the different pop- 
ulations as a function of the adopted magnitude cut. The re- 
sults are shown in Fig. 3 for both the standard initial mass func- 
tion — top panel — and the initial mass function of Adams & 
Laughlin (1996) — bottom panel. As it can be seen there, for 
the standard initial mass function the relative contribution of red 
dwarfs decreases with increasing magnitude cuts, while that of 
CO white dwarfs increases. Both contributions are equal for a 
magnitude cut of » 24"^^. Finally, the contribution of ONe white 
dwarfs remains roughly constant and only becomes appreciable 
for very large magnitude cuts. These trends can be attributed to 
the fact that red dwarfs are more numerous at bright magnitudes 
than white dwarfs, for which typical luminosities are of the or- 
der of log(L/L0) ^ -3.5. The situation is completely different 
when the log-normal initial mass function of Adams & Laughlin 
(1996) is used. As can be seen in the bottom panel of Fig. 3, 
the number of microlenses is practically dominated by the CO 
white dwarf contribution, while the contribution of red dwarfs 
and ONe white dwarfs is negligible. 

3.2. The microlensing event rate 

As previously pointed out, the contribution to the optical depth 
for a standard initial mass function is doubled when the red 
dwarf population is considered. We have also shown that for a 
standard initial mass function the contributions of red dwiirfs 
and CO white dwarfs are roughly the same. Then, it is natu- 
ral to ask ourselves if there are differences which can help us 
in discerning the contribution of one or another population us- 
ing the observational data of the MACHO experiment. To assess 
this we have analyzed the microlensing rate as a function of the 
event duration. The results of our simulations for a standard ini- 
tial mass function are shown in Fig. 4. Each of the panels is 
clearly labelled with the adopted magnitude cut and the popu- 
lation of microlenses. For the case of the left panels in Fig. 4 
we have adopted a magnitude cut of 25"^^, whilst for the right 
panels a magnitude cut of 30"^^ was adopted. In all the cases 
the simulated microlensing rate is shown using solid lines, while 
the observational data obtained by the MACHO team is shown 
using a dotted Une. All the distributions are normalized to unit 
area. As can be seen, the red dwarf and the white dwarf distribu- 
tions present some differences. For instance, for a magnitude cut 
of 25™''^ it turns out that although both the red dwarf population 
and the white dwarf populations show a peak located at nearly 
~ 20 days, the white dwarf population presents a wider distri- 
bution. In the case in which a magnitude cut of 30"^^ is adopted 
the differences are more pronounced and it is clear that the peak 
of the distribution for the CO white population moves to longer 



S. Torres et al.: The contribution of red dwarfs and white dwarfs to the halo dark matter 



7 



1 dr 
r. dt 



1 dr 
r„dt 



1 dr 
r.dt 




Fig. 4. Differential event rate normalized to unit area as a function of the Einstein crossing time for the populations under study and 
different magnitude cuts (soUd lines). Also represented as a dotted Une in each panel is the observational event distribution from 
Alcock et al. (2000). 
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Table 3. Compatibility, as obtained using the 7? statistical test, of the observed MACHO distribution of Einstein crossing times and 
those of the different simulated populations. 



durations (~ 70 days), whereas the peak of the red dwarf popu- 
lation does not move appreciably. 

Additionally, and in order to make more quantitative esti- 
mates, we have performed a statistical test of the compati- 
bility of the different populations under study with the observed 
data. The statistical test (Lucy 2000) represents and improve- 
ment over the standard statistical test and it is especially de- 
signed for meagre data sets. In Table 3 we show the 7? prob- 
ability of the different simulated populations being compatible 
with the distribution of Einstein times obtained by the MACHO 
experiment. It should be clarified that this probability represents 
an estimate of the degree to which the observed event rate distri- 
bution derives from a single population of stars. As can be seen 
in Table 3, the CO white dwarf population best matches the ob- 
servational data, given that its compatibility is as high as 0.90 
for the faintest magnitude cut. Moreover the compatibihty of 



this population with the observational data increases for fainter 
magnitude bins. In sharp contrast, the population of red dwarfs 
presents a decreasing trend as the magnitude cut increases and, 
additionally, the compatibility with the observational data is 0.70 
in the best of the cases. With regard to the ONe white dwarf 
population the compatibility presents and almost constant value 
around 0.70, independently of the magnitude cut. These results 
indicate that the CO white dwarf population can reproduce the 
observed distribution of microlensing event rates. Even more, 
they dominate the behavior of the whole population, as can be 
seen by taking a look at the last row of Table 3, in which we 
analyze the compatibility of the whole population of simulated 
stars. Therefore, even if the expected number of microlensing 
events obtained in our simulations is considerable smaller than 
the ~ 15 events claimed by the MACHO team, the event rate dis- 
tribution of the CO white dwarf population is in fair agreement 



8 



S. Torres et al.: The contribution of red dwarfs and white dwarfs to the halo dark matter 




-2.0 -1.0 0.0 1.0 



log(M/Mj 

Fig. 5. Halo dark matter fraction as a function of the mass lens. 
Plotted using a solid line are the 95% confidence level curve for 
the MACHO experiment and the 95% confidence level upper 
Umit for the EROS experiment. 

with the observed distribution. This result puts into question to 
which extent we know the characteristics of the halo white dwarf 
population and if it exists other possible ways capable of produc- 
ing a larger number of old white dwarfs in the stellar halo. 

3.3. The EROS experiment 

While the MACHO team claims for the identification up to 17 
observed events, the EROS collaboration have not found any 
microlensing event towards the LMC and one candidate event 
towards the SMC. Adopting a standard halo model and assum- 
ing TsMC = 1.4tlmc> the EROS results imply an optical depth 
To = 0.36 X 10"^ (Tisserand et al. 2007), which is four times 
smaller than that obtained by the MACHO team. We have also 
performed a set of simulations emulating the conditions of the 
EROS experiment with the same inputs previously described in 
Sect. 2. Although only minor differences should be expected in 
the analysis of the main results, it is clear as well that a joint 
study of both experiments using a controlled set of prescriptions 
represents a test of the robustness of our numerical procedure. 

In Table 4 we have summarized the results obtained in this 
second set of Monte Carlo simulations of microlenses towards 
the LMC for the EROS experiment. Our simulations show that, 
for the standard initial mass function, the expected optical depth 
could be 70% of the value found by the EROS team. The value 
obtained when only the white dwarf population was considered 
has been found to be 50% (Camacho et al. 2007). Thus, the sim- 
ulations presented here reproduce better the results of the EROS 
experiment. Obviously the red dwarf population is the responsi- 
ble of this increment. On the other hand, when a non-standard 
initial mass function is adopted, the results show only marginal 
differences with respect to those obtained for a white dwarf pop- 
ulation, given that in this case the role of red dwarfs is extremely 
limited. In summary, our results are in fair agreement with those 
obtained by the EROS experiment.and they seem to indicate that 
the microlensing optical depth obtained by the MACHO collab- 
oration is an overestimate. 



3.4. The dark matter density 

With the information obtained so far, we are in position to assess 
to which extent the mass range the stellar populations discussed 
here can contribute to the baryonic dark matter halo. Based on 
their ~ 15 microlensing events, the MACHO collaboration has 
derived an estimate of the halo fraction of dark matter / as 
well as the MACHO mass m using maximum-likelihood tech- 
niques. A similar analysis has been done by the EROS team, 
but with the significant difference that, in this case, no event 
has been reported for the LMC, which means that only an up- 
per limit on the halo mass fraction can be obtained. In order to 
compare the results of the MACHO and EROS collaborations 
with our Monte Carlo simulations we have adopted as our ref- 
erence model the isothermal sphere of core radius 5 kpc with a 
value of po = 0.0079 Mq pc"^ for the local dark matter den- 
sity and disregarding the contribution of the LMC halo. For 
this model we obtain that the optical depth towards the LMC 
is Tlmc = 5.1 X 10"^. The different estimates of the halo mass 
fraction / as a function of the mass are plotted in Fig. 5. We 
show as a solid line the curve of the MACHO 95% confidence 
level as taken from Alcock et al. (2000) as well as the EROS 
95% confidence level upper limit based on no observed events 
in the EROS-1 and EROS-2 data (Tisserand et al. 2007). Also 
represented are the individual contribution of each population 
under study as well as the whole population along with their 
corresponding 95% confidence level error bars. It is remarkable 
that the value obtained for the whole halo simulated population 
is in agreement within the 95% confidence level curves of both 
observational estimates. Thus, our results predict that the range 
of stellar masses within 0.08 and 10 M© provides / = 0.05 and 
an average mass of 0.41 1 Mq to the halo dark matter in agree- 
ment with the currently observational data. This result corrobo- 
rates our previous estimates about the limited contribution of the 
CO white dwarfs and the ONe white dwarfs (Garcia-Berro et al. 
2004; Camacho et al. 2007). 

4. Conclusions 

We have extended our previous studies about the contribution 
to the halo dark matter of the white dwarf population to the 
Galactic population of red dwarfs. Based on a series of Monte 
Carlo simulations which incorporate the most up-to-date evolu- 
tionary tracks for red dwarfs, CO white dwarfs and ONe white 
dwarfs and reliable models for our Galaxy and the LMC we have 
estimated the contribution of these objects to the microlensing 
optical depth towards the LMC and compared it with that ob- 
tained by the MACHO and EROS collaborations. In a first set 
of simulations we have found that the contribution of the red 
dwarf population practically doubles the contribution found so 
far for the white dwarf population. Our results indicate that the 
whole population of these stars can account at most for a ~ 0.3 
of the optical depth found by the MACHO team. This value im- 
plies that the contribution of the full range of masses between 
0.08 and 10 Mq represents a 5% of the halo dark matter with 
an average mass of 0.4 Mq. Despite that this result is in partial 
agreement with the 95% confidence level MACHO estimate for 
a standard isothermal sphere and no halo LMC contribution, the 
expected number of events obtained by our simulations (3 events 
at the 95% confidence level) is substantially below the 13 to 17 
observed MACHO events. These arguments reinforce the idea, 
previously pointed out by other studies, that the optical depth 
found by the MACHO team should be an overstimate, proba- 
bly due to contamination of self-lensing objects, variable stars 
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Table 4. Summary of the results obtained for the simulation of microlenses towards the LMC for the EROS experiment for an age 
of the halo of 14 Gyr, different model initial mass functions, and several magnitude cuts. 



and others. Moreover, we have assessed the compatibility be- 
tween the observed event rate distribution and the ones obtained 
for the different populations under study. Our results show that 
the CO white dwarf population can reproduce fairly well the ob- 
served event rate distribution although, as mentioned earUer, the 
expected number of events is considerable smaller. On the other 
hand, the negative results obtained by the EROS team towards 
the LMC are in agreement with our standard halo simulation. 
Finally, and for the sake of completeness, we have studied the 
effects of a log-normal biased initial mass function. In this case, 
the contribution of the red dwarf population is only marginal 
given that the production of low-mass stars is strongly inhib- 
ited. Accordingly, the total contribution to the microlensing op- 
tical depth is not different from that found in previous studies for 
the white dwarf contribution. 

Acknowledgements. Part of this work was supported by the MEC grants 
AYA05-08013-C03-01 and 02, by the European Union FEDER funds and by 
the AGAUR. 



References 

Adams, F.C., & LaughHn, G., 1996, ApJ, 468, 686 

Alcock, C, Allsman, R.A., Alves, D., Axelrod, T.S., Becker, A.C., Bennett, D.R, 
Cook, K.H., Freeman, K.C., Griest, K., Guern, J., Lehner, M.J., Mashall, 
S.L., Peterson, B.A., Pratt, M.R., Quinn, RJ., Rodgers, A.W., Stubbs, C.W., 
Sutherland, W., & Welch, D.L., 1997, ApJ, 486, 69 
Alcock, C, Allsman, R.A., Alves, D.R., Axelrod, T.S., Becker, A.C., Bennett, 
D.R, Cook, K.H., Dalai, N., Drake, A.J., Freeman, K.C., Geha, M., Griest, 
K., Lehner, M.J., Marshall, S.L., Minniti, D., Nelson, C.A., Peterson, 
B.A., Popowski, P, Pratt, M.R., Quinn, RJ., Stubbs, C.W., Sutherland, W., 
Tomaney, A.B., Vandehei, T., & Welch, D., 2000, ApJ, 542, 281 
Althaus, L.G., Garci'a-Berro, E., Isem, J., C6rsico, A.H., & Rohrmann, R.D., 

2007, A&A, 465, 249 
Baraffe, I., Chabrier, G., Allard, F, & Hauschildt, PH., 1998, A&A, 337, 403 
Becker, A.C. Rest, A. Stubbs, C, Miknaitis, G.A., Miceli, A., Covarrubias, 
R., Hawley, S.L., Aguilera, C, Smith, R.C., Suntzeff, N.B., Olsen, K., 
Prieto, J.L., Hiriart, R., Garg, A., Welch, D.L., Cook, K.H., Nikolaev, 
S., Clocchiatti, A., Minniti, D., Keller, S.C.& Schmidt, B.R, 2005, lAU 
Symposium, 225, 357 
Belokurov, V., Evans, N.W., Le Du, Y, 2003, MNRAS, 341, 1373 
Belokurov, V., Evans, N.W., Le Du, Y., 2004, MNRAS, 352, 233 
Binney, J., & Tremaine, H., 1987, Galactic Dynamics (Princeton: Princeton Univ. 
Press) 

Calchi Novati, S., de Luca, F, Jetzer, P, & Scarpetta, G., 2006, A&A, 459, 407 
Camacho, J., Torres, S., Isern, J., Althaus, L.G., & Garci'a-Berro, E., 2007, A&A, 
471,151 

Canal, R., Isern, J., & Ruiz-Lapuente, R, 1997, ApJ, 488, L35 
Chabrier, G., Segretain, L., & Mera, D., 1996, ApJ, 468, 21 
Chabrier, G., 2004, ApJ, 611, 315 
Chariot, S., & Silk J., 1995, ApJ, 445, 124 

Cutri, R.M., Skrutskie, M.F., van Dyk, S., Beichman, C.A., Carpenter, J.M., 
Chester, T., Cambresy, L., Evans, T., Fowler, J., Gizis, J., Howard, E., 
Huchra, J., Jarrett, T., Kopan, E.L., Kirkpatrick, J.D., Light, R.M., Marsh, 
K.A., McCallon, H., Schneider, S., Stiening, R., Sykes, M., Weinberg, 
M., Wheaton, W.A., Wheelock, S., & Zacarias, N., 2003, 2MASS All 
Sky Catalog of point sources, Univ. of Massachusetts and IPAC/Califomia 
Institute of Technology 



Dehnen, W. & Bmney, J., 1998, MNRAS, 298, 387 

Eisenstein, D.J., Liebert, J., Harris, H.C., Kleinman, S.J., Nitta, A., Silvestri, 
N., Anderson, S.A., Barentine, J.C., Brewington, H.J., Brinkmann, J., 
Harvanek, M., Krzesiliski, J., Neilsen, E.H., Long, D., Schneider, D.P., & 
Snedden, S.A., 2006, ApJS, 167, 40 
Finley, D.S., Koester, D., & Basri, G., 1997, ApJ, 488, 375 
Flynn, C, Holopamen, J., & Hohnberg, J., 2003, MNRAS, 339, 817 
Garcfa-Ben-o, E., Tortes, S., Isem, J., & Burkert, A., 2004, A&A, 418, 53 
Gates, E.I., & Gyuk, G., 2001, ApJ, 547, 786 

Goldman, B., Afonso, A., Alard, Ch., Albert, J.-N., Amadon, A., Andersen, J., 
Ansari, R., Aubourg, E., Bareyre, P., Bauer, F., BeauUeu, J.-Ph., Blanc, G., 
Bouquet, A., Char, S., Chariot, X., Couchot, F., Coutures, Ch., Derue, F., 
Ferlet, R., Fouqufi, R, GUcenstein, J.-F, Gould, A., Graff, D., Gros, M., 
Halssinski, J., Hamadache, C, Hamilton, J.-Ch., Hardin, D., Kat, J. de, 
Kim, A., Lasserre, Th., Le Guillou, L., Lesquoj, E., Loup, C, MagnevUle, 
Ch., Mansoux, B., Marquette, J.-B., Maiuice, E., Maury, A., Milsztajn, A., 
Moniez, M., Palanque-Delabrouille, N., Perdereau, O., Prfivot, L., Regnault, 
N., Rich, J., Sprro, M., Tisserand, P., Vidal-Madjar, A., Vigroux, L., & 
Zylberajch, S., 2002, A&A, 389, 69 

Griest, K., 1991, ApJ, 366, 412 

Gyuk, G., Dalai, N., & Griest, K., 2000, ApJ, 535, 90 

Hambly, N.C., Irwhi, M.J., & MacGilUvray, H.T., 2001, MNRAS, 326, 1295 

Hansen, B.M.S., 1998, Nature, 394, 860 

Hansen, B.M.S., & Liebert, J., 2003, ARA&A, 41,465 

Holopamen, J., Flynn, C, Knebe, A., GUI, S.P., & Gibson, B.K., 2006, 

MNRAS368, 1209 
Ibata, R.A., Richer, H.B., GiUiland, R.L., & Scott, D., 1999, ApJ, 524, L95 
Iben, I., & LaughUn, G., 1989, ApJ, 341, 312 

Isem, J., Garcla-Berto, E., Hemanz, M. Mochkovitch, R., & Torres, S., 1998, 

ApJ, 503, 239 
James, F, 1990, Comput. Phys. Commun., 60, 329 

Kallivayalil, N., van der Marel, R.R, Alcock, C, Axelrod, T., Cook, K.H., Drake, 

A.J., & Geha, M., 2006, ApJ, 638, 772 
Kerins, E., & Evans N.W., 1998, ApJ503, 75 

KUic, M., Mendez, R.A., Von Hippel, T., & Winget, D.E., 2005, ApJ, 633, 1126 
Klypm, A., Zhao, H., & Somerville, R., 2002, ApJ, 573, 597 
Knox, R.A., Hawkms, M.R.S., & Hambly, N.C., 1999, MNRAS, 306, 736 
Lasserre, T. Afonso. C, Albert, J.N., Andersen, J., Ansari, R., Aubourg, E., 
Bareyre, P., Bauer, F., Beauheu, J.P., Blanc, G., Bouquet, A., Char, S., 
Chariot, X., Couchot, F, Coutures, C, Deme, F, Ferlet, R., GUcenstein, 
J.F, Goldman, B., Gould, A., Graff, D., Gros, M., Haissinski, J., Hamilton, 
J.C., Hardin, D., de Kat, J., Kim, A., Lesquoy, E., Loup, C, Magneville, 
C, Mansoux, B., Marquette, J.B., Maurice, E., Milsztajn, A., Moniez, M., 
Palanque-Delabrouille, N., Perdereau, O., Pr^vot, L., Regnault, N., Rich, J., 
Spho, M., Vidal-Madjar, A., Vigroux, L., & Zylberajch, S., 2001, A&A, 
355, L39 

Liebert, J., Bergeron, R, & Holberg, J., 2005, ApJS, 156, 47 

Liebert, J., Dahn, C.C., & Monet, D.G., 1989, m "White Dwarfs", Ed. G. Wegner 

(Berhn: Springer), 15 
Lucy, L., 2000, MNRAS, 318, 92 
Majewski, S.R., & Siegel, M.H., 2002, ApJ, 569, 432 

Mollerach, S., & Roulet, E., 2002, Gravitational Lensing and Microlensing 
(Singapore: World Scientific) 

Muraki, Y, Sumi, T., Abe, F, Bond, I., Carter, B., Dodd, R., Fujimoto, M., 
Heamshaw, J., Honda, M., Jugaku, J., Kabe, S., Kato, Y, Kobayashi, M., 
Koribalski, B., Kilmartin, P., Masuda, K., Matsubara, Y, Nakamura, T., 
Noda, S., Pennycook, G., Rattenbury, N., Reid, M., Saito, T., Sato, H., 
Sato, S., Sekiguchi, M., Sullivan, D., Takeuti, M., Watase Y, Yanagisawa, 
T., Yock, P. & Yoshizawa, M., 1999, Progress of Theoretical Physics 
Supplement, 133, 233 

Navarto, J.F, Frenck, C.S., & White, S.D.M., 1997, ApJ, 490, 493 



10 



S. Torres et al.: The contribution of red dwarfs and white dwarfs to the halo dark matter 



Nelson, C.A., Cook, K.H., Axelrod, T.S., Mould, J.R., & Alcock, C, 2002, ApJ, 

573, 644 

Oppenheimer, B.R., Hambly, N.C., Digby, A.P., Hodgkin, S.T., & Saumon, D., 

Science, 292, 698 
Paczynski, B., 1986, ApJ, 304, 1 
Rahvar, S., 2005, MNRAS, 356, 1127 
Reyle, C, Robin, A.C., & Creze, M., 2001, A&A, L53 
Sahu, K., 1994, Nature, 370, 275 

Salaris, M., Garcfa-Berro, E., Hemanz, M., Isem, J. & Saumon, D., 2000, ApJ, 

544, 1036 

Scalo, J., 1998, in "The Stellar Initial Mass Function", Eds.: G. Gilmore & D. 
Howell (San Francisco: PAS? Conference Series), Vol. 142, 201 

Schneider, P., Kochanek, C.S., & Wambsganss, J., 2004, in "Saas-Pee lectures 
on Gravitational Lensing" 

Silk, J., 2007, in "The invisible universe: dark matter and dark energy", Proc. of 
the Third Aegean Summer School (Berlin: Springer Verlag), Lecture Notes 
in Physics, Vol. 720, 101 

Tisserand, P., Le Guillou, L., Afonso, C, Albert, J. N., Andersen, J., Ansari, 
R., Aubourg, E., Bareyre, P., Beaulieu, J. P., Chariot, X., Coutures, C, 
Ferlet, R., Fouque, P., Glicenstein, J. F., Goldman, B., Gould, A., Gralf, D., 
Gros, M., Haissinski, J., Hamadache, C, de Kat, J., Lasserre, T., Lesquoy, 
E., Loup, C, Magneville, C, Marquette, J. B., Maurice, E., Maury, A., 
Milsztajn, A., Moniez, M., Palanque-Delabrouille, N., Perdereau, O., Rahal, 
Y. R., Rich, J., Spiro, M., Vidal-Madjar, A., Vigroux, L., & Zylberajch, S., 
2007, A&A, 469, 387 

Torres, S., Garci'a-Berro, E., & Lsern, J., 1998, ApJ, 508, L71 

Torres, S., Garci'a-Berro, E., Burkert, A., & lsern, J., 2002, MNRAS, 336, 971 

Udalski, A., Szymanski, M., Kaluzny, J., Kubiak, M., Mateo, M., & Krzeminski, 
W., 1994, Acta Astron., 44, 1 

Vennes, S., Smith, R.J., Boyle, J., Croom, S.M., Kawka, A., Shanks, T, Miller, 
L., & Loaring, N., 2002, MNRAS, 335, 673 

Vidrih, S., Bramich, D.M., Hewett, PC. Evans, N.W., Gilmore, G., Hodgkin, 
S., Smith, M., Wyrzykowski, L., Belokurov, V., Fellhauer, M., Irwin, M.J., 
McMahon, R.G., Zucker, D., Munn, J.A., Lin, H., Mrknaitis, G., Harris, 
H.C., Lupton, R.H., & and Schneider, D.R, 2007, MNRAS, 382, 515 

Zhao, H.S., 1998, MNRAS, 294, 139 



